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ABSTRACT 


In  this  paper  presentj^the  results  of  a limited  number  of  ex- 
periments with  linear  fractional  problems.  Six  solution  procedures 
were  tested  and  the  results  are  expressed  in  number  of  simplex-llke 


pivots  required  to  solve  a sample  of  twenty  problems  randomly  generated. 


ap 


Two  main  approaches  emerge  from  Che  literature  to  solve  the 
linear  fractional  problem; 

V-  max  { f (x)“n(x)/d{x)  : x CF}  (P) 

where  n(x)“c^+c  x , d(x)“d^+dx,  F»  { x e R*'  : Ax“b  x^O}, 

c and  d are  real  numbers,  c and  d are  real  n-vcctors,  A is  an  mxn 
o o 

real  matrix  and  b is  a real  m-vector.  We  assume  in  this  note  that  F 

Is  compact  and  Chat  min  { d(x)  : x c F }>  0. 


Charnes  and  Cooper  (A)  transform  problem  (P)  into  the  linear  program: 

V-  max  { c^t+cy  : Ay-bt»0,  d^t+dy-1,  and  t,y^  0}. 

This  approach  has  been  extended  to  the  nonlinear  version  of  (P)  by  Bradley 
and  Frey  {3]  and  Sohaible  [7).  The  second  approach  solves  a sequence  of 
linear  problems  or  at  least  one  pivot  step  of  each  linear  program  over 
the  original  feasible  set  by  updating  the  objective  function.  Algorithms 
in  this  category  are  related  to  ideas  first  presented  by  Isbell  and  Marlow  (5) 
and  Martos  [61.  Similar  algorithms  have  been  proposed  by  several  other 
authors.  The  interested  reader  is  referred  to  Che  excellent  bibliography 
collected  by  I.M.  Stancu-Mlnasian  [8].  Methods  in  the  second  approach  pro- 
pose to  solve  (P)  through  a sequence  of  linear  programs: 

r(x*^)-  max  { r(x*^,x)-n(x)-f (x^'ldCx)  : x e F } k-0,1,2...  (LPj^) 
where  x**  is  a given  feasible  point  and  x*'  for  k > 1 is  defined  in  Isbell  and 

Marlow's  procedure  as  being  the  optimal  solution  to  and  as  Che  first 

k-1 

feasible  basis  in  which  c(x  ,x)  >0  in  Martos'  procedure.  Both 

algorithms  tervlnate  at  iteration  k^  for  which  r(x  °)-0.  In  this  case 

It  is  worth  noting  that  Wagner  and  Yuan  [9]  related  the  two 

optimal. 

main  approaches  by  showing  that  Marcos'  algorithm  is  equivalent  to  Charnes  and 


Cooper's  method  in  Che  sense  that  both  algorithms  lead  to  an  identical  se- 
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quence  of  pivot  tri)!  operations.  Bltran  atiJ  Magnanti  [1]  have  extended  the 
connection  between  these  approaches  by  relating  th«'in  to  generalized  program- 
• Bing.  No  theoretical  or  empirical  evidence  has  been  given,  in  the  past,  in- 

dicating which  of  the  several  existing  algorithms  is  to  be  preferred. 

In  this  note  we  present  the  results.  In  number  of  simplcx-llke  pivots, 
of  twenty  problems  of  type  (P),  randomly  generated,  solved  by  the  following 
six  algorithms  (each  problem  when  solved  by  each  of  the  six  procedures  was 
I started  with  the  same  basic  feasible  solution): 

* 

A)  Maximize  n(x)  over  the  feasible  set  obtaining  the  optimal  solution  x . 

o * , 

NoxC«  apply  Isbell  and  Marlow  s algorithm  with  x ”x 

* 

B)  Minimize  d(x)  over  the  feasible  set  obtaining  the  optimal  solution  x . 

o * 

Next,  apply  Isbell  and  Marlow  s algorithm  with  x ^x  . 

t 

I C)  Maximize  g(x)"(c-(cd/dd)dlx  over  F obtaining  the  optimal  solution  x . 

o * 

[ Next*  Apply  Isbell  and  Marlow  s algorithm  with  x "x  (Bitran  and  Novaes 

I (21  suggested  the  objective  function  g(x)). 

[ 

I D)  Isbell  and  Marlow's  algorithm. 

L 

I E)  Harto  s ' algorithm. 

[ 

I F)  The  author  considered  relevant  to  compare  these  algorithms  with 

I the  number  of  pivots  necessary  to  solve  the  linear  programs: 

1 max  { n(x)-vd(x)  : x eF}  (LP) 

I where  for  each  of  the  twenty  problems  (P) , v Is  chosen  as  Its  optimal  value. 

^ The  optimal  value  of  (LP)  Is  zero  and  any  solution  to  (LP)  Is  optimal  In  the 

fractional  program  (P)  (ll)).  Note  that  (LP)  corresponds  to  (LPj^)  with 


’‘‘‘‘^optimal. 
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, The  characteristics  of  the  data  of  the  twenty  randomly  generated 

problems  are  the  following: 

n“40,  m“20,  the  absolute  value  of  each  • the  (l,J)th  element  of 

each  matrix  A was  randomly  generated  in  the  Interval  (0,10].  The  density  of 

negative  elements  being  207..  Each  component  b^  l>l,2.,,,m  of  each  right  hand 

n 

side  b was  defined  as  7 a.,/2.  The  objective  function  coefficients  c ,c  ,d 

IJ  o j o 

djj“l,2..n  were  generated  in  the  Intervals  [-1000  Cj  ^ 1000;  0<d^,  dj  < 1], 

[Kc  , c.  < 1000;  1<  d , d,  < 2]  or  [-1000  < c , c < -1;  1 < d , d < 2] . The 
“o’j-  -..o’j-'  -o’J-  ’-o’J-' 

reason  for  choosing  such  Intervals  was  to  obtain  five  problems  with  an  angle  0 be- 

cd 

tween  the  gradients  of  the  numerator  and  denominator,  l.e.,  cos  9 » pclflid^ 

In  each  of  the  four  Intervals  [O,]!),  [it,  tt)  , [ _^ ) , [ 3ti^  , 2 it  ) in  an  at- 

4 4 2 2 4'  4 

tempt  to  Identify  a correlation  between  the  algorithms  tested  and  the  geometry 
of  linear  fractional  programs.  The  geometric  properties  of  problem  (P)  are  con- 
sequences of  the  following  facts. 

I)  The  hyperplanes  n(x)-Ld(x).r0  contain  for  each  L both  the  sets  {xCr”: 
f(x)»L  } and  CE"  {x  G r”:  n(x)“0  and  d(x)“0}.  The  set  CE  Is  called  the 
center  of  the  problem  because  as  L varies  the  hyperplanes  rotate  about  it 
giving  a "star"  centered  at  CE([2)^ 

II)  The  objective  function  f(x)  Is  pseudo-concave  and  quaslconvex  on  the 

set  {x  C r”:  d(x)  >0}  , l.e.,  f(y)  >f(x)  If  and  only  If  7 f (x) (y-x) >0 . 

2 

In  R the  geometry  of  (P)([21)  suggests  that  procedure  (C)  would  perform 
better  than  (A)  and  (B)  for  high  and  low  values  of  9 (9  G (0,n]).  Table  1 
shows  the  results  obtained.  For  the  first  and  last  five  problems  a total  of  178 
pivots  was  necessary  with  procedure  (C)  while  233  and  363  pivots  were  requlrea 
with  procedures  (A)  and  (B)  respectively.  The  corresponding  standard  deviations 


r 
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being  3.70,  6-01  and  7.90.  For  the  twenty  probUos  selected  Hartos* 
algo  rlthm  performed  better  than  the  preceding  four  and  in  some  cases  required 
fewer  pivots  than  procedure  (F) . Algorithms  (C)  and  (0)  were  practically 
equivalent  and  were  followed  by  (A),  while  (B)  performed  poorly.  The  com- 
puter code  used  to  solve  the  twenty  problems  by  the  six  algorithms  was  an 
adaptation  of  Burroughs'  conmercial  code  TEMPO. 


I 
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